Rare predicted loss-of-function variants of type I IFN immunity genes are associated with critical COVID-19 (preprint)

Background We previously reported inborn errors of TLR3- and TLR7-dependent type I interferon (IFN) immunity in 1-5% of unvaccinated patients with life-threatening COVID-19, and autoantibodies against type I IFN in another 15-20% of cases. Methods We report here a genome-wide rare variant burden association analysis in 3,269 unvaccinated patients with life-threatening COVID-19 (1,301 previously reported and 1,968 new patients), and 1,373 unvaccinated SARS-CoV-2-infected individuals without pneumonia. A quarter of the patients tested had antibodies against type I IFN (234 of 928) and were excluded from the analysis. Results No gene reached genome-wide significance. Under a recessive model, the most significant gene with at-risk variants was TLR7, with an OR of 27.68 (95%CI:1.5-528.7, P=1.1x10-4), in analyses restricted to biochemically loss-of-function (bLOF) variants. We replicated the enrichment in rare predicted LOF (pLOF) variants at 13 influenza susceptibility loci involved in TLR3-dependent type I IFN immunity (OR=3.70 [95%CI:1.3-8.2], P=2.1x10-4). Adding the recently reported TYK2 COVID-19 locus strengthened this enrichment, particularly under a recessive model (OR=19.65 [95%CI:2.1-2635.4]; P=3.4x10-3). When these 14 loci and TLR7 were considered, all individuals hemizygous (n=20) or homozygous (n=5) for pLOF or bLOF variants were patients (OR=39.19 [95%CI:5.2-5037.0], P=4.7x10-7), who also showed an enrichment in heterozygous variants (OR=2.36 [95%CI:1.0-5.9], P=0.02). Finally, the patients 13 with pLOF or bLOF variants at these 15 loci were significantly younger (mean age [SD]=43.3 [20.3] years) than the other patients (56.0 [17.3] years; P=1.68x10-5). Conclusions Rare variants of TLR3- and TLR7-dependent type I IFN immunity genes can underlie lifethreatening COVID-19, particularly with recessive inheritance, in patients under 60 years old.

Genome-wide detection of human variants that disrupt intronic branchpoints (preprint)

ABSTRACT Pre-mRNA splicing is initiated with the recognition of a single-nucleotide intronic branchpoint (BP) within a BP motif by spliceosome elements. Fifty-six rare variants in 44 human genes have been reported to alter splicing and cause disease by disrupting BP. However, until now, no computational approach has been available to efficiently detect such variants in next-generation sequencing (NGS) data. We established a comprehensive human genome-wide BP database by integrating existing BP data, and by generating new BP data from RNA-seq of lariat debranching enzyme DBR1-mutated patients and from machine-learning predictions. We in-depth characterize multiple features of BP in major and minor introns, and find that BP and BP-2 (two-nucleotides upstream of BP) positions exhibit a lower rate of variation in human populations and higher evolutionary conservation than the intronic background, whilst being comparable to the exonic background. We develop BPHunter as a genome-wide computational approach to systematically and efficiently detect intronic variants that may disrupt BP recognition in NGS data. BPHunter retrospectively identifies 48 of the 56 known pathogenic BP mutations in which we summarize a strategy for prioritizing BP mutation candidates, and the remaining 8 all create AG dinucleotides between BP and acceptor site which is probably the reason for mis-splicing. We demonstrate the utility of BPHunter prospectively by using it to identify a novel germline heterozygous BP variant of STAT2 in a patient with critical COVID-19 pneumonia, and a novel somatic intronic 59-nucleotide deletion of ITPKB in a lymphoma patient, both of which we validate experimentally. BPHunter is publicly available from https://hgidsoft.rockefeller.edu/BPHunter and https://github.com/casanova-lab/BPHunter .

The risk of COVID-19 death is much greater and age-dependent with type I IFN autoantibodies (preprint)

SARS-CoV-2 infection fatality rate (IFR) doubles with every five years of age from childhood onward. Circulating autoantibodies neutralizing IFN-α, IFN-ω, and/or IFN-β are found in ~20% of deceased patients across age groups. In the general population, they are found in ~1% of individuals aged 20-70 years and in >4% of those >70 years old. With a sample of 1,261 deceased patients and 34,159 uninfected individuals, we estimated both IFR and relative risk of death (RRD) across age groups for individuals carrying autoantibodies neutralizing type I IFNs, relative to non-carriers. For autoantibodies neutralizing IFN-α2 or IFN-ω, the RRD was 17.0[95% CI:11.7-24.7] for individuals under 70 years old and 5.8[4.5-7.4] for individuals aged 70 and over, whereas, for autoantibodies neutralizing both molecules, the RRD was 188.3[44.8-774.4] and 7.2[5.0-10.3], respectively. IFRs increased with age, from 0.17%[0.12-0.31] for individuals <40 years old to 26.7%[20.3-35.2] for those ≥80 years old for autoantibodies neutralizing IFN-α2 or IFN-ω, and from 0.84%[0.31-8.28] to 40.5%[27.82-61.20] for the same two age groups, for autoantibodies neutralizing both molecules. Autoantibodies against type I IFNs increase IFRs, and are associated with high RRDs, particularly those neutralizing both IFN-α2 and -ω. Remarkably, IFR increases with age, whereas RRD decreases with age. Autoimmunity to type I IFNs appears to be second only to age among common predictors of COVID-19 death.